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Abstract 
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I. INTRODUCTION 



Interactions of atoms with strong electric and magnetic fields, in particular crossed fields, 
have long been the focal point of various research efforts [???????]. The 
conflicting symmetries inherent in the electronic motion in atomic and external fields make 
for a rich variety of phenomena. In this light they are not only interesting from a theoretical 
standpoint, such as the hydrogen atom, which constitutes a paradigm for a chaotic system 
[?????]: Their study also had a major impact on a number of other fields, ranging 
from semiconductor physics to astrophysics. 

However, for a long time theoretical investigations did not go beyond the stage of treating 
the nucleus as 'frozen', the so-called infimte-mass approximation. The model underlying this 
picture is the concept that the atomic motion could be separated into that of the center of 
mass and the one relative to it. While this is perfectly justified for translationally invariant 
systems — going along with the conservation of the total momentum — it ceases to be true 
in the presence of magnetic fields. Nonetheless, a weaker version of this procedure can be 
established in the case of (neutral) atoms in homogeneous fields [????]. In this pseudo- 
separation^ a conserved quantity termed pseudo-momentum — the total momentum plus a 
field-dependent compensation — allows to obtain a Hamiltonian that depends on the relative 
coordinates only. However, that relative motion strongly depends on the center of mass via 
the pseudo-momentum. This coupling gives rise to some new and interesting effects, such 
as chaotic diffusion of the center of mass [? ? ], or dynamical self-ionization [? ? ] for the 
case of ions. 

Yet one of the most prominent among these finite-mass effects is the existence of giant 
dipole states, where the electrons can be decentered from the nucleus by many 10, 000 atomic 
units [? ? ]. The systematic base for this was the gauge-independent generalized potential 
first derived for hydrogen by Dippel et al. [? ] and extended to arbitrary atoms by 
Schmelcher [? ]. It has been applied extensively to study giant dipole states in the two- 
body case; e.g., hydrogen [? ] and positronium [? ? ], where the large inter-particle 
distance prevents annihilation for up to several years. By contrast, the extension to the 
multi- electron case has attracted only limited attention. There have been indications for 
two-electron quasistable giant-dipole states based on a local semiclassical analysis [? ]. 
The aim of this paper is to both provide a local analysis of general A^-electron system and 
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investigate the stability of the full quantum-mechanical system. In order to first analyze 
the A^-electrons case analytically, we look for stationary points of the so-called generalized 
potential of the effective relative motion (the 'giant-dipole configuration') so as to carry 
out a normal-mode analysis about these points. Equipped with this insight into the local 
behavior, we study the exact system N = 2 numerically using wave-packet propagation. 
Here we resort to the mctdh method (Multi-Configuration Time-Dependent Hartree) [? ? 

This article is organized as follows. In Sec. Ull the theoretical framework of the pseudo- 
separation is reviewed and applied to derive the generalized potential of the relative motion. 
Sec. mil then deals with the stationary points of the generalized potential as a base for the 
analysis of giant dipole resonances. In the subsequent section, the normal-mode analysis 
for electrons is carried out. To this end, the local equations of motion are derived 
()IV Ajl ■ whose eigenmodes and -vectors are computed numerically (jIV Bjl as the solutions of 
a quadratic eigenvalue problem. Sec. El contains a wave-packet dynamical study of the two- 
electron system based on the mctdh method. Results on the stability, spectral properties 
and lower bounds for the lifetimes of the resonances are presented. 

II. THE iV-ELECTRON ATOM IN CROSSED ELECTRIC AND MAGNETIC 
FIELDS 

The Hamiltonian of an atom, consisting of N electrons (of mass m) and a nucleus (mass 
Mo) interacting via the Coulomb potential V, reads in the laboratory frame: 

TV N 

^ = ^ E - + ^ + ^^Mr'o)? - e ^ E ■ r: + ZeE ■ r'^ + V. 

i=l i=l 

Here an electron (index i < N) has the electric charge e = — |e|; the atomic number Z equals 
the number of electrons N in the neutral case treated here. We consider static electric and 
magnetic fields, which are accounted for by the electro-static potential 0(r) = — E ■ r and 
the vector potential A(r), whose gauge is not fixed here. 

In the field-free case, the translation invariance of the system — going along with the 
conservation of the total momentum — guarantees a complete separation of the center-of- 
mass and the relative motion. The vector potential now breaks that symmetry. For the 
special case of homogeneous fields, though, it was shown that a so-called pseudo-separation 
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of the center-of-mass motion is possible . The key is that even though the total momentum 
is not conserved, the total pseudo-momentum K := '^^=0^^ defined in terms of the one- 
particle pseudo-momenta 

kj := mvj — er'j x B 

— and likewise for the nucleus — is constant by construction [???]. It can be thought of as 
the kinetic center-of-mass (CM) momentum plus a magnetic-field-dependent compensation, 
which trivially reduces to the total momentum in the case B = 0. Moreover, its components 
{Ka} mutually commute for a neutral system, making them commensurable constants of 
motion. This is crucial in performing a gauge-independent pseudo-separation of the CM 
motion [? ? ], which consists of the following steps: 

1. Transformation to CM and coordinates relative to the nucleus, {RcMj^i}, yielding 
HcM and Kcm- 

2. Choose the total wave function to be a simultaneous eigenstate of the pseudomomen- 
tum with eigenvalue K: (Kcm — K)\I/k = 0. The result is of the form 

^k(Rca./, {r,}) = UiRcM, {r,})0({r,}), (1) 

where is an arbitrary function of the relative coordinates, and U denotes a unitary 
operator. 

3. Finally reduce the full Schrodinger equation to an effective one for the relative coor- 
dinates via 

{n - E)(t) = 0, with UTiK-^ = Hcm- (2) 

It has to be emphasized that gauge independence is an essential quality of the pseudo- 
separation scheme outlined here. A major drawback inherent in any fixed-gauge approach is 
that one cannot identify gauge-independent terms in the effective Hamiltonian making 
the analysis less systematic. 

In order to derive this effective Hamiltonian in a gauge-invariant way, we decompose the 
vector potential according to 

A(r) = -B X r + VA(r). 
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The transform U and hence the effective Hamiltonian Ti. — T -\-V is then determined only 
up to some A-dependent function /({rj}) of the relative coordinates which vanishes by 
construction for VA = [? 1: 



2m 9 



N 

m 



i=i 



2 

+ eV,/| + 



i(E[p.-^Bxr, + eV,/]) 



(3) 




^ ^ 2i(^^""^^>^"^ l -e>:E.r, + l^, (4) 

with q = Are(Mo - Nm)/M. 

T represents the (gauge-dependent) kinetic energy of the relative motion; it depends 
on the gauge via /. As opposed to that, the generalized potential V is manifestly gauge- 
independent and is therefore identified as an effective potential for the electronic motion at 
a given value of the total pseudo-momentum. Beyond the usual Coulomb potential and the 
Stark terms due to the external electric field, we have an expression corresponding to the 
CM-kinetic energy, 

— (^K-eBxgr.j ^ — + -[b .^.^r,^ (B) 

The first term, the analogue of the field-free CM energy, is merely a constant. The second 
expression is a diamagnetic term, which acts as a confining harmonic potential on the elec- 
trons' center of mass, j;jYl!j=i'^i-> with a frequency ijOqm = NeQ/M perpendicular to the 
magnetic field. Most notably, the center of mass — unlike in the field-free case — actually cou- 
ples to the electrons (i.e., to their center of mass) via the motional electric field E^,/ := ■ 
For crossed fields, it is therefore inviting to combine the Stark terms arising from both 
the external and the motional electric field and account for them by defining the effective 
pseudomomentum 



K' := K + M——- = K - Mu, 



E X B 

u being the classical drift velocity. The Hamiltonian then reads [? ? 



n{K, E) = n{K', 0) + K • u - ^Mu^. 
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We have thus reduced our problem of crossed fields to a purely magnetic one. In what 
follows, we will always deal with TC{K', 0) for simplicity. 

The generalized potential derived above gives rise to many effects, among the most promi- 
nent ones being the existence of giant dipole states (GDS). The concern of the following two 
sections is to derive stationary points of that potential which accomodate these states, and 
to classify their normal modes. 

III. STATIONARY CONFIGURATIONS OF THE GENERALIZED POTENTIAL 

With the gauge-independent generalized potential at hand, we can now start our search 
for decentered stationary points, which we will identify as candidates for giant dipole reso- 
nances. To facilitate the search for stationarities, we will first choose convenient coordinates 
(Subsection IIII Ajl which present a suitable extension of the case = 2 [? ]. We will then 
derive and solve the stationarity conditions by an educated guess in IIII Bl 

A. The electronic center of mass frame 

As we know from the studies on the one-electron case [? ? ? ], it is the CM-kinetic energy 
^ that is responsible for the decentered states. It was pointed out in Eq. (jH) already that 
this term depends on the combination R = ;^X]j=i^j' which we refer to as the electronic 
center of mass (ECM). It is thus natural to introduce the ECM as a new coordinate. In 
addition, we will choose A^ — 1 vectors relative to the ECM, i.e. to decompose = R + Sj, 
where the {sj} only account for all but one (say, the A^th) electron's relative position: 

(ri,...,rAr) t-^ (R;Si,...,SAr_i) 

1=1 

1 ^ 

Si:=ri-— ^r^- {i = 1, . . . , N - 1), 
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Introducing the conjugate momenta P = —i-^, t^j 
becomes (with / = for simplicity) 

1 



-i^, the Hamiltonian Ti = T + V 



r 



2fiN 2 



BxR)2 + ^y L-ly 

^ 2m ^ I N ^ 

i<N \ k<N 



TTfc - -B X Sj + 



i<N 



i<N 



V 



a=l a<b<N 



(6) 
(7) 



where we introduced the shorthands 

Sat = FAf - R = - ^ Sj 



i=l 



ra = kal = |R + Sal 
Sab = |Sa - Sfel {a <b < N) 

as well as the reduced mass yU^^ = {Nm)~^ + Mq"^. It is reassuring that the kinetic energy 

can be shown to have the intuitive form 

1 . 1 ^ 
r = -/i7vR^ + -^ms^. 

a=l 

B. Existence of stationary points 

Given the advantageous form of the Hamiltonian (jMIj), we are now ready to look for 
stationary points of the potential. If we define Q = (R, Si, . . . , SAr_i)^ G M^^, then we seek 
some Qo such that 

dV 

This stationarity condition leads to the rather intricate set of equations, 
A^eB 



dV 



M 



X (K' - A^eB X R) + 



(8) 



i<N 
1 1 



i<N 



N 



'i ' N/ 'i 'TVJ 



(9) 



El 1 V ■\ S^. V ■\ / I 1 V ■\ 1 

\ k<N ^ik ^iN k<N '^ik f.<N \ j<N j'^ 

\ kj^i I k^i 
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(Throughout this subsection, the subscript in Qq is omitted for mnemotechnical reasons.) 

Equations ^ and Q constitute a 3A^-dimensional non-linear system of equations that 
has little hope of being solved analytically in complete generality. However, one can try to 
find special solutions by an educated guess. For instance, one immediately sees that Eq. ^ 
simplifies considerably if r, = for all i. At that stage, all electrons are distributed on 
the shell of a sphere centered around the nucleus. It is tempting to go one step further and 
require R ± Sj so as to 'decouple' the two types of coordinates. This motivates our set of 
assumptions I: 

= r^v =: r, R ■ Si = Vz < N. (10) 

Before summarizing the immediate consequences, let us specify a basis rj, ^} such that 
the magnetic field points along the ^ direction: 

C = B/S 

rj = K'/K' 

^ = ry X C/\V X CI- 

(For convenience, we assume Z(K', B) = 90° in addition to the premise of crossed fields. 
This may be achieved by preparing a state with K ~ so that only the external electric 
field is relevant.) In this basis, write the extremal ECM vector in terms of R-^ = X, R-ry = 
Y,R-C = Z. 

If (|TTHl holds, then the extremal condition (jH)), = 0? is equivalent to 

1. Y = Z = 0, that is R = {X, 0, 0)'^, where 

2. X satisfies P,(X) = + ^-^X + NZe^^ = 0. 

3. Moreover, there is a set of angles (pi such that Sj = s (0, cos0j, sin0j)^ Vi < A^. 

By (1.), if there is a solution, then the electronic center of mass will be aligned along the 
$, axis, that is the decentering takes place anti-parallel to the electric field. Based on this, 
statement (2.) gives the (implicit) condition for the existence of solutions X . Finally, (3.) 
refiects a constraint on all relative coordinates s^. they have the same norm s and — due 
to their orthogonality to R — can be placed arbitrarily on a circle about R. For symmetry 
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reasons, we demand that all relative coordinates be distributed uniformly over the circle, 
which amounts to the assumption II: 

27r 



N' 



(11) 



This procedure determines the circular configuration only up to a global rotation by an angle 
A G [0,27r). 

The above guess allows us to fix the ratio of |X| and s, or r and s, respectively, so as to 
meet Eq. ©: 



AN 



r = s. 



E.<^(smf)-i 



r = (l-aj/) 2 |X| =: aAr|X|. 



(12) 
(13) 



Pictorially, is related to the angle 5n between R and via sin^Tv = cnj^^. To be explicit, 
we have computed the parameters for some relevant N below: 





aN = i 




2 


2 


^ =1.15470054... 


3 


VS = 1.7320508... 


= 1.22474487... 


4 


1.6107688... 


1.27558847... 


5 


1.5372471... 


1.31666491... 



The connection between r and X now also enables us to turn the implicit equation 
Pr{X) = above into a cubic equation 

K' M 



(14) 



Eq. (HH) has two distinct solutions Xi/2{K; B) provided that K' > K^r = ^ 



Xx = 7;^(2cos( " ' : " • 1-1), A = 1,2 (15) 



cos 6* 



3NB 



K' 



1. 



In other words: for any B there is a critical value Kcr for K' below which |^(X^, sf') = 
has complex solutions. In an experiment, one therefore has to control the value of K' = |K + 
ME X B/i?^|. Averbukh et al. [? ] suggested how this could be achieved experimentally. 
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Figure 1: The extremal component ^1/2 as a function of K/Kcr for different field strengths i3/a.u. 
The critical point X{Kcr) corresponds to the vertex where Xi (lower graph) and X2 (upper graph) 
meet. 

The underlying idea is to prepare an atom with low K (requiring a slow CM motion), and 
use the electric field E to control the critical parameter. This way an initially prepared 
Coulomb-Rydberg state is transformed into a decentered state localized in the extremum. 

To gain a better understanding of the two solutions, their behavior is illustrated in Figure 
m For fixed S, both lines meet at the vertex = Xi{Kcr) = X2{Kcr) = —Kcr/NB. From 
there on, the lower branch Xi{K) tends to —00 monotonically (which we identify with a 
proper decentering), whereas the upper curve X2{K) goes to zero (i.e., the decentering gets 
lost). For stronger magnetic fields, in turn, the 'decentering region' around Xc shrinks. We 
remark that the decentering is virtually insensitive to the number of electrons, N . 

Let us now wrap up our results. Every stationary electronic vector possesses the orthog- 
onal decomposition 



This means that all electrons have a common position along ^ — which we refer to as the 
decentering — while they are confined to a highly symmetric circular configuration perpen- 
dicular to ^. This circular configuration in turn is determined only up to an overall rotation 
A. A generic setup is shown in Figure [21 

The extremal solution exists only if the effective pseudomomentum K' (controlled via 
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Figure 2: A giant-dipole configuration for iV = 4 electrons {B = W~'^, K = 2Kcr). The ECM 
R = {X, 0, 0) is decentered, with the relative vectors Sj confined to a circle as indicated. Also 
shown is the electronic vector relative to the nucleus rj = R + Sj. A typical example of the 
cyclotron and the motion due to a Coulomb mode is indicated (see Subsec. IIV B Hi . 

the electric field) exceeds some critical value. The more it does so, the more distinct the 
decentering of the ECM will be. 

IV. NORMAL-MODE ANALYSIS OF iV-ELECTRON GIANT DIPOLE STATES 

As argued in the previous section, v^e expect the stationary configurations to be candi- 
dates for giant dipole resonances With certain lifetimes. Preceding a numerical study, we 
first seek to obtain some insight into the local stability of the system, i.e., their behavior in 
the vicinity of the extrema. To that end, v^^e will carry out a normal-mode analysis about 
these points. More specifically, we shall consider the harmonic system 

V{Q) := V{Q^) + ]^{Q - Qofw{Q ~ Qo), 

with the Hessian matrix w = 9^V((5o) G M^^^'^^. According to Ehrenfest's theorem, in this 
case the expectation values {Q)(t) obey the corresponding classical equations of motion. 
Their solution leads us to a quadratic eigenvalue problem, which will be solved by lineariza- 
tion so as to obtain the eigenmodes and -vectors, giving an insight into both the dynamics 
and the spectrum of the problem. 
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A. Equations of motion 



Before stating the equations of motion, we are to deduce the form of the Hessian matrix 
of the generalized potential V. Using the obvious notation 



N-l 







aaaR aRas^ 
a^v a^v 



as^aR asias,- 



WRR WRj 
WiR Wij 



a lengthy calculation reveals that the symmetric matrix w has the the 3x3 block structure 
{NeBf 



Wrr 



M 



(16) 



WRi = -Ne^— (2S {(si - Sat) O R} + O Sj - Sjv O Sat) 



,2_ 

r-5 



-Ne^— (R R(l + 5ij) + 2S {(si% + sjv) (g) R} + (0 Si% + sjv (0 sjv) + 



AT 



3 

kN 



a=l 



3e- 



N 



(Sj - S„) ® {Sj - Sa) 



(sjv - Sfe) (8) (siv - Sfc) 



+ 



+ 



a=l 



fe<Af 



(si - stv) O (si - Sn) 



(1 - Sij) 



^Nk 

[Si - Sj) (g) {Si - Sj) 



)■ 



where /-*- = diag(l, 1,0), S denotes the symmetrizer: S{A} — ^(^4 + ^4^), and {<S>] represents 
the tensor product. (All coordinates here of course refer to the extremal positions.) 

Despite this rather involved form, the equations of motion themselves look fairly simple: 



/xjvR = qRxB -wr - {Q - Qo) 
eSj X B - 



mSi 



(wi - -^^w^ ■ {Q - Qo). 



(17) 



Here a mixed notation was employed for clarity: Wa = {waR, Wai, ■ ■ ■ ,Wa,N-i) £ M^^^-^ is 
the a-th block row of w. It is inviting to interpret these equations as the motion of effective 
particles with masses hm, tti and charges g, e, respectively, in a magnetic field, coupled via 
linear forces. 

To solve the equations of motion, let us first write the above system in a more compact 
way in terms of the displacements v{t) :— Q{t) — Qo, 



V = uu ■ V + A ■ V. 



(18) 
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A 




3Nx3N 



Here the antisymmetric cyclotron matrix 00 and the harmonic matrix A, respectively, have 
been introduced: 

eB eB 
— e © ■ ■ ■ © — 

m m 

{wi -jfj^k^k) /m 

where {e)ij := Eij^ is connected with the Levi-Civita tensor and [©] stands for the direct sum. 
The solution of this 3iV-dimensional second-order system of differential equations (ODE) is 
given by the span 

6N 

<t) = Y.^v,e^^')c, {v, G C^^; 7p, G C), (19) 
p=i 

fulfilling the quadratic eigenvalue equation 

(7,'/-7p^-^)^P = 0. (20) 

In this respect, our stability analysis amounts to finding the complex eigenvalues 7p =: 
T p + iVtp (whose imaginary parts are frequencies of a vibration about a stable point, and 
whose real part corresponds to an instability), and the eigenvectors Vp. 

The above quadratic eigenvalue problem is solved via linearization, which is entirely 
equivalent to reducing a second-order ODE to a first-order system. In this fashion one 
obtains a standard linear eigenvalue problem 

Au = 7M, (21) 

with the linearization matrix 

' / 



A 



A u 

Before we present the results, let us state two general properties of the solutions that can 
be obtained ahead of a computation: 

1. For any eigenpair (7,f) of (|2n|l such that Im7 7^ , there is another pair (7*,^*). 

2. Let {uj)y := v^'uv denote the Rayleigh quotient for an arbitrary matrix u. Assuming 
f to be a normalized eigenvector, then for any mode the eigenvalue is given by 



^ 2 

where in general only one of the roots is in the spectrum. 
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Figure 3: Eigenmodes {7p}p for = 2 electrons as a function of the electric field E = BK/M 
over the range K/Kcr S [1,10] (all modes are imaginary). The three plots refer to the magnetic 
fields B = 10"^, 10"^, 10~^ (from left to right). Each top horizontal line represents two almost 
degenerate cyclotron modes. Below, the two Coulomb modes fall off quickly and intersect the CM 
mode (the nearly horizontal line about four orders below the cyclotron modes). 

Firstly, (1.) implies that only one out of two modes is relevant; hence it is widely legitimate 
to treat the 6A^ modes as effectively 3A^. (2.) tells us where in the complex plane we can 
expect the eigenmodes to lie. If the problem is elliptic, that is 

(^).<-(^)' v., (22) 
then the solution is imaginary 7 = if2 (since uj is antisymmetric). 

B. Results and discussion 

1. General Classification 

Let us first have a look at the behavior of the modes {^p{K)}p for different at fixed B 
(Figures 01 SJ . Even though the patterns become increasingly rich and involved with higher 
A^, one can see a distinction between the modes regarding their behavior as a function of K 
(and B), their order of magnitude and, more generally, their location in the complex plane. 
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Figure 4: Eigenmodes 7 for the cases = 3 (a), = 4 (b), = 5 (c) and = 10 (d), as a 
function of the electric field E = BK/M, plotted over the range K/Kcr G [1,10]. The magnetic 
field is = 10~^. Imaginary parts appear as solid lines, while points (+) are used for the real 
parts. The individual modes are assigned just as for Fig. 01 with an extra class of very small decay 
modes, which fall off twice as fast as the Coulomb lines. 

An analysis of the associated eigenvectors and the criterion for ellipticity, (f^ . v^^ill also 
allovi^ us to illuminate the modes' origins. This suggests the following classification of the 
3A^ eigenmodes: 

• N so-called cyclotron modes corresponding to the cyclotron motion of the effective 
particles. Their values are exclusively on the scale of the cyclotron frequencies, 7 = 
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Figure 5: Visualization of the zero mode: plot of the generalized potential V for = 2 close to 
the extremum, cut along the yz coordinates of s [K/Kcr = 2; B = 10~^). There is a circle that 
continuously connects all minima, referred to as 'stationary circle'. 

•n -lel-B 

m 

• Coulomb modes corresponding to the inter-Coulombic motion, with a main contri- 
bution from the harmonic matrix A. They fall off quickly with K, since the decentering 
Xq increases, and so the Coulomb interaction becomes less effective. 

• 1 CM mode, roughly reflecting the cyclotron motion of the center of mass , O, ~ 
(stemming mostly from the CM-kinetic energy) 

• 1 zero mode (7^ = 0) that stems from the rotational invariance of the saddle point 
with respect to the circular configuration (see Eq. fTT| 

• N — 2 modes that will be referred to as decay modes in recognition of the fact that 
they are predominantly real. They are neither directly related to the cyclotron motion 
nor to the spectrum of A, and their slope is twice as steep as that of the Coulomb 
modes. 

The cyclotron modes roughly pertain to gyrations perpendicular to B, while the Coulomb 
motion is predominantly parallel to B. The CM mode is primarily a gyrational mode — 
reflecting the nucleus motion — , but it also exhibits Coulomb-type contributions. The zero 
mode — to be analyzed in more detail below — corresponds to the rotation along the stationary 
circle identiflable in Fig. El Lastly, the 'decay modes' reveal no striking underlying structure. 
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There are a few subtleties that go beyond the categorization suggested above. A closer 
look at, for instance, the case of = 10 electrons reveals certain interactions among the 
different modes. Their principal causes are crossings between the center-of-mass mode and 
the decay modes (resulting in some striking deformations of the usual line pattern), and 
equally avoided crossings of the CM mode with at least some of the Coulomb modes. This 
may be taken as a hint at the different symmetry relations among the Coulombic and the 
decay modes. 

S. Influence of the rotational invariance A 

By construction, the stationary character was not affected by a common rotation of the 
relative coordinates Sq by an angle A on the circle: 

Sa ^ -DaSq Va, (23) 

4>a ^ 0a + A 

where Da G 1 © S0(2) denotes that rotation about the ^ axis. Nonetheless, the Hessian 
matrix 9^VL = w is affected — and thus the dynamics differs. To see this, note that the 
block matrices transform like 

wrr ^ Wrr {N > 2) 

Waj ^ DAWajD^K 

The spectra of both the Coulombic Hessian {d'^V)o and the magnetic contributions (includ- 
ing u) are separately invariant under these rotations. Still, their conflicting symmetries 
guarantee that the dynamics will vary with A. 

In order to see this dependency, a semi-log plot of the modes {7p(A)} on the maximal 
interval [0, 27r) is recorded for different A^ (Fig. As before, B = 10""^ and K = 2Kcr are 
fixed. The plots display an obvious symmetry, as rotating by $7v = ^n/N gives an indis- 
tinguishable setup. The cyclotronic modes as well as the affiliated CM mode are essentially 
A-independent. By contrast, the Coulomb modes and the decay modes show a pronounced 
periodic change. In addition, one can see avoided crossings in some cases between the CM 
mode and the lowest Coulomb mode. 

For the case N = 2 (Figure EK), which was found to be locally stable in the vertical 
configuration ([? ], see also lIVBl) . there is some exceptional behavior whenever A is close to 
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(b) Case of three electrons 

Figure 6: Upon rotating the circular configuration by an angle A, the spectrum is changed in 
part {B = 10-"^, K = 2Kcr). The cyclotronic modes (658 GHz) and CM mode (0.17GHz) remain 
constant. The Coulomb modes and, for = 3, the decay mode (~ lO^^GHz) show some periodic 
modulation, (for symbol coding, see Fig. 0J 
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Figure 7: Spectrum for = 4 depending on the rotation angle A {B = W , K = 2Kcr); see text. 



k ■ Tc. The lower Coulomb mode has an avoided crossing with the almost constant CM mode 
before it experiences a singularity on the logarithmic scale, along the way turning real. In 
this respect, the local stability we found for = 2 is applicable only outside the singular 
horizontal configuration A = 0. 

While the case of three electrons (A^ = 3: Fig. IHb) reveals no striking effects, the spectrum 
becomes much richer when we add another electron (N = 4: Fig. [Zj). The four Coulomb 
modes as well as two more decay modes reveal an intriguing behavior. There is an avoided 
crossing between the CM mode and the lowest Coulomb mode, which turns into an unstable 
real mode in an intermediary region. Also, one can get a sense of the interplay between the 
two decay modes. 



V. WAVE-PACKET DYNAMICAL STUDY OF THE TWO-ELECTRON CASE 

The extensive analysis of the normal modes of A^-electron giant dipole states in the 
preceding section provided evidence that the two-electron case is locally stable, apart from 
the singular horizontal configuration. We now turn to a numerical study of that system. 

It cannot be over-emphasized that a six-dimensional resonance study is at the frontier 
of what is currently possible and requires a careful choice of the computational approach. 
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This applies especially in view of the fact that the system is governed by dramatically dif- 
ferent time scales. Therefore we adopted the Multi-Configuration Time-Dependent Hartree 
(mctdh) method [????], which is known for its outstanding efficiency in high- 
dimensional applications. To be self-contained, a concise introduction to this tool is pre- 
sented in Subsection IV Al Its basic idea is to solve the time-dependent Schrodinger equation 
in a small but time-dependent basis related to Hartree products. It has to be stressed that 
this method is designed for distinguishable particles. Applying it to a fermionic system like 
ours finds its sole justification in the fact that the spatial separation between the electrons 
is so large that they are virtually distinguishable. 

After some remarks on the scales of the system in Subsec. IV B| the method will be applied 
in Subsec. IV ("I The focus is on stability aspects, i.e., the propagation of wave packets, but 
we also investigate the spectra as well as the infiuence of the rotational invariance discussed 
in the previous section. 

A. Computational method: MCTDH 

The principal goal of mctdh is to solve the time-dependent Schrodinger equation as an 
initial- value problem: 

= Hm 

(24) 

*(Q,0) = vI/o(g) 

The standard method (see [? ] and references therein) now approaches the problem by ex- 
panding the solution in terms of time-independent primitive basis functions. Unfortunately, 
it exhibits a drastic exponential scaling with respect to the number of degrees of freedom /. 
Rather than using a large time-independent basis, the Multi- Configuration Time-Dependent 
Hartree method (mctdh) employs a time-dependent set, which is more fiexible and thus 
smaller. The ansatz for the wave function now reads 

vi/(g,t) = 5^...^A,,..,,(t)ny,(:)(g.,t)^ J2 Ajmj{Q,t), 

ji=l jf=l K=l \J\<\n\ 

using a convenient multi-index notation for the configurations, J = (ji...j/), where / 
denotes the number of degrees of freedom. The (unknown) single-particle functions (pf^ 
are in turn represented in a primitive basis. Note that both the coefficients Aj and the 
Hartree products $j carry a time-dependence; hence the uniqueness of solutions can only 
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be ensured by adding certain constraints on the single-particle functions [? ? ]. Using the 
Dirac-Frenkel variational principle [? ], one can derive equations of motion for both $j 
[? ? ]. Integrating these ODE systems allows us to obtain the time evolution of the system, 
^('i*^) = ^)) which lies at the heart of wave-packet propagation. 

The Heidelberg mctdh package [? ], which we used, incorporates a few extensions to 
this basic concept: 

• Product representation of the potential: in order to circumvent multi-dimensional in- 
tegration in computing the matrix elements, MCTDH makes the requirement 

/ 



enforcing that the numerical Hamiltonian be written as a sum of products of one- 
particle operators (separable form). The non-separable part of the Hamiltonian thus 
has to be fitted to product form ahead of a computation. 

• Mode combination: in practice, one combines several (say d — 2,3) degrees of freedom 
to a d-particle rather than a one-particle function. This approach alleviates the bad 
numerical scaling for high-dimensional systems. 

• Relaxation [? ]: MCTDH provides a way to not only propagate a wave packet, but 
also to obtain the lowest eigenstates of the (discretized) system. The underlying idea 
is to propagate some wave function by the non-unitary time-evolution operator e~^'^ 
{propagation in imaginary time.) In practice, one relies on a more sophisticated scheme 
termed improved relaxation. Here {'^\H — El"^) is minimized with respect to both the 
coefficients Aj and the configurations $j. The equations of motion thus obtained are 
then solved iteratively by first solving for Aj{t) (by diagonalization of {{^j\H^k)) 
with fixed $j) and then propagating $j in imaginary time over a short period. The 
cycle will then be repeated. 

• Spectrum: one can also extract information on the spectrum by computing the auto- 
correlation function c{t) :— {^q \ e~'^^*^Q). In the case of a purely discrete spectrum, 
its Fourier transform reads 




(25) 




(26) 



J 
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that is, all eigen-energies give a peak-type contribution, weighted according to the 
overlap cj of \&o with the J-th eigenstate. 

B. Application of MCTDH 

Let us point out how the computational method can be applied to our problem. There are 
essentially two highly conflicting types of motion and corresponding scales: perpendicular 
to the magnetic field, we have the magnetic length 



employing atomic units. To be specific, for a laboratory field strength B = 10~^, the length 
scale will be on the order of 10^ a.u. As for the parallel (Coulomb) motion — expected to 
take place roughly in a harmonic potential — we can use the usual harmonic-oscillator lengths 
Zq = 1/^/JEoz, for some 'typical' frequency Uz- These are roughly on the order of 10^ a.u. 
(Note that the higher K — or E — becomes, the smaller will be the Coulomb frequencies, 
hence the wave packet will be spread far more over the grid.) 
Analogously, the anticipated time scales are 

Tb = — ^ 1.5 ps 

Ub 

for the cyclotronic motion. By contrast, the characteristic times parallel to B are strongly 
i^T-dependent, and they vary between about 1, 500ps— 15, OOOps. For comparison, we mention 
that the associate time scale for the CM mode is asymptotically Tcm — 5, OOOps. 

For strong fields, the relevant motion will thus certainly be determined by the magnetic 
field — at least near the stationary points we are interested in. This would suggest to in- 
troduce cylindrical coordinates about the extremum, carrying out our study in the basis 
set for the Landau orbitals [? ] Rnmip)^^"^''^ perpendicular to B, and harmonic-oscillator 
functions 4>n:,{z) parallel to the field. Since the former set was not available in the mctdh 
package, we resorted to Cartesian coordinates for the displacements from the extremum Qq 
as a pragmatic solution: 

s s — So = (x, y, 2;), R 1-^ R — Ro = (X, Y, Z). 

The Hamiltonian (EHU) for our case, centered about the extremum Qo by the generalized 
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translation operators expi(k/j • Rq + ■ Sq), reads: 

V = ^ ((2e5)2 [X^ + Y^] + 2X ■ {Xo{2eBf - 2eBK')) +V0 + V1 + V2 + C 

where ujr = \q\B/fi, Ug = \e\B/m. C = ^(K.' — 2eB x Rq)^ only shifts the energy by a 
constant and will be ignored hereafter, and vq = <5|Xo|, Sq = r^/a, 0o = + A denote the 
extremal configuration. Moreover, Vq stands for the electron-electron repulsion, whereas V1/2 
denotes the attractive potentials between the nucleus and electron #1 and #2, respectively. 
Both T and the 'CM-kinetic energy' are already in product form (f^ . whereas the Coulomb 
terms have to be fitted. 



C. Results and discussion 



In our calculations, we were chiefly interested in two questions: the stability of the giant 
dipole states, and their spectral properties. The flrst aspect was studied using propagation 
of wave packets initially localized at the extremum. Emphasis was placed on initial states 
close to the assumed 'ground states'. Moreover, following the evolution of wave packets 
with an initial displacement from the extremum, we tested the robustness of the resonances. 
For comparison, the eigenvectors of the Hamiltonian will be considered using the relaxation 
method. 

The parameter set was restricted to the following values. B = 10"^ was fixed, because 
changing the magnetic field strength essentially only affects the overall scales via the critical 
pseudomomentum and is not expected to make for completely new behavior. The pseudo- 
momenta are K/Kcr G {1.1, 2.0, 10.0} (we will drop the prime from now on), which account 
for the cases of a very shallow outer well, the medium range and the very deep outer well 
(the high-i^ regime). 

For simplicity, the influence of the rotational invariance will be scrutinized a posteriori, 
using the characteristic values 0o = f ^ ^5 ^- These refer to the vertical configuration 
(A = — I) first studied in [? ], a 'diagonal' one A = —j, and the horizontal configuration 
A = 0, which we expect to have an instability according to the analysis in IIV B 2l 
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The propagation times were chosen in the regime of 5 — 10 x 10'^ ps. This includes many 
periods of the Coulomb modes and several 10, 000 periods of the rapid cyclotron motion. 

1. Stability analysis 

The approach used in this subsection is to study the propagation of a 6D harmonic- 
oscillator wave packet centered within the well. The oscillator parameters are adapted so as 
to match the borderline ground states, i.e., the cyclotron frequencies ur, for XY{xy), and 
a typical Coulomb frequency uz, uj^ parallel to the magnetic field. To test the robustness 
of the potential resonance, in some cases displacements of {Z)q = {z)o = 2000 a.u. are also 
applied to the initial wave function. (This value is arbitrary but moderate compared to the 
total decentering.) Also, relaxations (starting from the original, undisplaced wave packets 
above) are carried out to see how the lowest-energy configuration found with this algorithm 
differs from our initial guess. The analysis mostly focuses on the reduced densities (ID or 
2D) obtained by integrating out all but the K-th degree(s) of freedom. Where appropriate, 
the wave packets' centers and widths {qK)(t), AqK{t) are shown. 

a. The case K = l.lKcr-' Very close to the critical point, K^r, the observed motion 
displays an instability in some degrees of freedom. 

Both X and Y revealed weak oscillations of two types. As these are generic, they shall be 
inspected in more detail on the example of the X degree of freedom. Three representative 
snapshots (Figure [HK) were taken to illustrate this effect. The nature of the oscillations 
can be explored more easily in Figure [Htb). To begin with, the center carries out a tiny 
oscillation about X ~ —20 a.u. with a period of about 6,400ps (a snapshot was taken at 
its first turning point at 3, 200 ps). This very slow oscillation superimposes what we will call 
the shape oscillations, materializing in the fiuctuation of the width AX(t). That oscillation 
is about twice as fast, and it reaches its maximum width for the first time at l,600ps 
(Fig. IHt^)- The relative stability in these two directions is readily explained in terms of 
the generalized potential. Apart from a stabilization by the B field, they experience an 
additional confinement via the 'CM-kinetic energy' (0). However, it is evident that the 
packet broadens over time, thus slowly delocalizing (Fig. |Hb)- 

Even though the relative motion perpendicular to B, = {x,y), is also gyrationally 
stabilized, there is no confining term for it in the generalized potential as for the ECM. In 
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(a) (b) 
Figure 8: The motion in X in the case K/Kcr = 1.1 

(a) Snapshots of the motion px , reflecting oscihations of both center and shape of the wave packet 

(b) The two oscillations as documented in {X){t), AX(t) 





Tc{K) 


\Xo{K)\ 


1.1 


.1 - .4 ns 


1.4- 10^ a.u. 


2 


.4 - 1 ns 


2.9 • 10^ a.u. 


10 


5 — 15 ns 


1.5 • lO^a.u. 



Table I: Summary of time scale of the Coulomb motion (Tc) and length scale (|-Yo|) for different 
values of pseudo-momentum K 

fact, a look at the generalized potential (Figure IHb) reveals that Qq is a saddle point. The y 
direction is the zero mode, while the cut through the x direction refers to a maximum. In this 
light, we cannot expect the system to be arbitrarily stable in these two degrees of freedom. 
This is a general fact, but for K = l.lKcr it is very pronounced. The reason being, near the 
extremum, the Coulomb interaction (the source of these instabilities) becomes negligible for 
higher K. A line of reasoning closely connected uses the time scale at which the instability 
will have an impact. As the Coulomb modes Q{K) increase dramatically as K ^ their 
characteristic time T{K) = 27t/Q{K) becomes much shorter (cf. Table HI). 

These considerations shed light on Figure [TTl For instance, px soon develops two small 
side maxima. They become more pronounced as time goes on and slowly spread towards the 
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(a) (b) 

Figure 9: Cuts through the generahzed potential near the extremum,V(X, Y, 0; 0) (a) and 
V(O;x,y,0) (b), exemplified for K = l.lKcr- The (XY) potential is almost harmonic. The y 
direction virtually coincides with the zero mode, whereas the potential exhibits an instability with 
respect to x. 



Generalized Potential in meV Generalized Potential in meV 




(a) Cut through the Z direction (b) Cut through the z direction 

Figure 10: Cuts through the generalized potential parallel to the magnetic field: (a) V(0, 0, Z; 0) is 
roughly harmonic, (b) V(0; 0, 0, z) is strongly flattened on the right-hand region of the minimum. 
(Note that the extremum is shifted to the origin.) 

boundary of the grid. The zero mode (y) spreads out much faster — that is, to a significant 
extent after 2,000 ps already. 

The reduced density for Z (Fig. HSJ broadens and continually leaks toward the grid 
ends. Although in the Z direction, the particle lives in a fairly harmonic well (Fig. ITHll . it 
is not only lifted in energy to way above the bottom of the well, but also affected by the 
instability in x via coupling. Turning to the last degree of freedom, z, we find an overall 
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(c) (d) 

Figure 11: Snapshots of the wave packet's motion in the direction of x (a), y (b) and z (c) for 
K = l.lKcr- The plots display the reduced densities px, Py and pz at different times. While the 
unstable mode x leaks slowly, the zero mode y spreads much more quickly. The wave packet in z 
is only marginally shifted over time. In (d), the different oscillations in z are illustrated. 

stable behavior. Note that there is an overlay of a fast but tiny oscillation of the packet's 
center and a slower and less regular one visible in Figure fTlT d) as the envelope of {z)(t). 

b. The case K = 2Kcr: For K being twice its critical value, the vertical configuration 
is virtually stable on a scale of T ~ 10"^ ps. The degrees of freedom X and Y, as well 
as the parallel motion Z and z, turn out to be indeed perfectly stable. Again, there are 
two different types of tiny oscillations, displayed in Figure [131 which are fingerprints of the 
coupling with the Coulombic motion (indicated by plotting {qK)it) over some periods Tz, 



27 




Figure 12: The motion of the Z mode in the case K/Kcr = l-l- Three snapshots of the reduced 
density pz are plotted. 




time / ps time / ps 

(a) (b) 

Figure 13: The expectation values and uncertainties over selected time periods for the coordinates 
(a) Z and (b) z {K/Kcr = 2). 

for instance). 

As opposed to the previous paragraph K/Kcr = 1-1, the relative motion in x conveys a 
fairly stable if erratic impression. Still, the wave packet in y spreads as it did before (Fig. 
mi). But keep in mind that it occurs at a much greater time scale than in Fig. ITlT b): for 
example, after 5,000 ps the broadening has only just begun to become significant. 

A closer look at the response of the system upon displacing Z and z by 2, 000 a.u. unveils 
that for the excited degrees themselves, Z and z, the wave packet is simply reflected between 
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Figure 14: The relative coordinate s±_ in the case K/Ka- = 2: Characteristic snapshots of the 
reduced density py. 

two positions ±(-2)0 with minor (Z) or more pronounced {z) deformations and smearing-out 
due to competing modes, cf. Figure fT5lf a-b). The most interesting question may be the 
effect on the {xy) modes. In fact, x is rendered slightly unstable by the excitation of the 
parallel motion. The snapshots taken at 20,000 ps and 50,000 ps, respectively, document 
how the wave packet slowly but inevitably starts leaking, Fig. fTH c-d) . The lifetime of the 
y mode is not drastically changed altogether. 

We also investigated the supposed ground state obtained via improved relaxation. It 
is understood that at best the exact system displays a resonance state, so the quest for 
eigenstates is meaningful only insomuch as those of the discretized Hamiltonian may be 
interpreted as localized states which really correspond to resonances. Our results turned out 
to be very sensitive to both the primitive and the single-particle basis sizes. Nonetheless, for 
a high if manageable accuracy, the pattern that emerged was the following. Apart from some 
smaller deviations from the initial state (Landau-orbital / harmonic-oscillator product), a 
drastic alteration takes place in s± = {x,y) (Fig. fT6|l . The wave packet is split into two 
humps that seem to be driven outwards until the basis size is exhausted, while smeared out 
along the zero mode y. 

This behavior raises the question how that squares with the potential picture, given that 
the local instability characterized by = Wxx/'fn is much weaker in magnitude than a;^. In 
order to clarify these issues, we constructed a 2D model Hamiltonian via the natural inclusion 
Hxy '■= H |pi=o;z=o- The study of this toy model, which is numerically far more amenable, 
revealed that it is capable of recovering many of the key features of the full system. It 
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Figure 15: Snapshots of the wave-packet upon an initial displacement hy Z = z = 2,000a.u. for 

K/Kcr = 2. Plot (a) illustrates the strong deformations in Z in the course of the oscillation; (b) 
shows how the the z packet is distorted after several thousand pico-seconds. The non-stable modes 
x,y are affected via coupling, see (c) and (d). 

therefore provided a valuable tool in detecting the key to this mechanism — the paramagnetic 
term cub^z- To begin with, the instability of the generalized potential in x tempts the packet 
to split up and be driven 'downhill'. However, if there were no paramagnetic term or if 
the same instability occurred isotropically, this would have no discernible effect. However, 
owing to the anisotropy of V, £ is no longer conserved: 

£ = -s X VV. 

It is energetically favorable for the particle to go to ever lower (i^) < while being expelled 
from the unstable line x — 0, thus introducing a high correlation between x and y. 
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Figure 16: A generic relaxation state {K = 2Kcr, see text) - plotted is the 2D-reduced density for 
the plane perpendicular to B: the wave packet is driven away from the x axis with increasing basis 
size. 

Note that this effect essentially amounts to the fact that the giant-dipole resonance state is 
far from the energetically lowest configuration. That does not pose an immediate constraint 
on the system's lifetime under propagation, since our initial Landau/harmonic-oscillator 
state has a virtually zero overlap with the above 'relaxation state'. Hence in a propagation, 
this 'bimodal' relaxation state is not accessible, which explains why the system is so much 
more stable than its relaxation suggests. Physically speaking, the cyclotron gyration field 
stabilizes the motion near the stationary configuration for fixed energy and thus inhibits 
a delocalized state at a lower energy. - This is approved by a 2D numerically exact time 
evolution (Fig. ITTjl . which reveals no breakup for many nanoseconds, until the wave packet 
is slightly rotated in the xy plane and split up roughly in the y direction — and not in x as 
in the relaxation case! Note that it acquires a small but non-zero value of {i^) ~ —1, in 
contrast to ~ —60 in Fig. HU this also reflects the higher level of stability. It is only 
after way more than 100 ns that the packet hits the grid's boundary. These outcomes seem 
to confirm that for moderate values of K, our system is quasi-stable on a time scale of up 
to T ~ 100 ns. 

c. The case K = lOKcr: For this value, which corresponds to E = BK/M = 207^, 
we find the system to be practically stable on the time scale we considered, T = 100 ns. 

The initial wave packet experiences only tiny deformations in most degrees of freedom, 
which are not discussed here. However, there are oscillations for z, which are more pro- 
nounced but still marginal compared to the spatial extension of the decentered state; these 
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Figure 17: Wave packet of the 2D model (xy) after T = 160 ns (K = 2Kcr^ see text) 
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Figure 18: The time evolution of {z),Az (a) and Ay (b) for K/Kcr = 10. Both degrees are 
altogether stable, but oscillate due to coupling. A marginal broadening can be seen in y. 

are sketched in Figure [THa). While its center shows a shift to positive values of its elon- 
gation is on the order of 250 a.u., and the width varies by a factor of two, all at a period 
of about Tz. In contrast with x and any other degree, the width in the y direction shows a 
steady if minute broadening (Fig. fHflh). This might be taken for a sign that eventually, the 
system is bound to decay. Then again, this takes 300, 000 ps to move by one atomic unit 
and is beyond the time scale regarded here. 

After having come a long way to find the system practically stable for high enough values 
of the pseudo-momentum, we now explore the robustness with respect to displacements. 
Skipping the detailed behavior, we find that the overall stability was not affected. Due to 
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Figure 19: The circular configuration s^'') = so(0, cos (po, sin (^q)'^ for N = 2, depending on the angle 
00- The arrows refer to s^*^) in the yz-plane for different values of 0o. The X axis points into the 
page. 

coupling to other modes, the excitation of the parallel motion imprints some oscillations 
both of the packet's center and its width on the perpendicular degrees of freedom. 

Lastly, applying a relaxation revealed similar phenomena as in the previous paragraph. 
Regardless of the specific shape of the relaxation state, it can be reasoned that the insta- 
bilities of course still underlie the system, but that the lifetimes for sufficiently large K are 
too large to be observed in our propagation. 

2. Influence of the rotational freedom 

To complete the discussion, let us now investigate the giant dipole states for different 
settings of the extremal angle 0o of the relative coordinate Sq on the circle. So far, the focus 
was on 00 = 7r/2 only (or, equivalently, A = — f), which we referred to as the 'vertical 
configuration'. In the normal-mode analysis, however, we found an interesting frequency 
pattern as we varied this extremal angle. It was argued that each 0o corresponds to a 
different extremal point, i.e., for each value of {B,K) there is really a one-parameter class 
of extremal points Q{4>o) with different stability and spectral properties. 

As examples, we thus investigate both the supposedly unstable 'horizontal configuration' 
00 = TT (A = 0) and a 'diagonal configuration' 0o = = as visualized in Figure 

Uni This was done for the special cases of K/Kcr = 2, 10. 
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Figure 20: The horizontal configuration (/)q = tt {K/Kcr = 10). The wave packet switches between 
two positions (snapshots of reduced density for 11,500 ps and 20,000 ps). 

To sum up our findings, the horizontal configuration indeed adds an instability, which 
is discernible even for very high K, if less distinct. The diagonal configuration was partly 
unstable on a timescale comparable to that of the vertical case, T ~ 10^ ps. A rotation of 
the extremum obviously rotates the zero mode, too, which is why the detailed dynamics is 
a different one for the relative degrees of freedom. For instance, in case that (po = tc, the 
zero mode points along the z direction, which would otherwise live in a harmonic potential. 
As a consequence, even for K/Kcr = 10 the z-packet (prepared with a non-zero frequency 

= 10~^a.u.) smoothly oscillates between its initial and a much more smeared-out state 
according to Fig. 1201 

3. Spectrum 

We finally look into some spectral properties of the GDS for the different values of K/Kcr- 
The spectrum calculated from a propagation for K/Kcr = 1-1 (Fig. 1211) looks somewhat 
fuzzy, but the main peak is still distinct enough to vaguely resemble a resonance. The 
situation is less ambiguous in the case K/Kcr = 2 , as the 'decay' occurs on a much greater 
time scale. The spectrum in Fig. l22l illustrates the resonance character of the system. Our 
initial wave packet gives a peak at ~ 0.281 meV, whereas the 'displaced state' with 
{z)o = {Z)q = 2, 000 a.u. produces a rather interesting excitation spectrum (same figure), 
revealing overlaps with many eigenstates. The equidistant spacing of the peaks might be 
interpreted as a signature of harmonicity in both excitations {Z, z). Finally, in the most 
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Figure 21: Fourier transform of the auto-correlation function (the 'spectrum') for K/Kcr = l-l- 
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Figure 22: Spectrum for K/Kcr = 2. The dotted peak again belongs to the initial state. For the 
excited (or displaced) initial state, the spectrum shows a sequence of peaks (here almost equidis- 
tant). 

stable case K/ Kcr = 10, the spectrum consists of a sharp dominant peak at about 2.2422 eV, 
plus a very tiny one 0.0017 eV above. This somehow quantifies our empirical observation 
that K/Kcr — 10 is quasi-stable. 

VI. CONCLUSION AND OUTLOOK 

We have studied the internal motion of A'"-electron atoms in crossed magnetic and electric 
fields utilizing the generalized potential. For electric fields above a critical value, a strongly 
decentered potential well forms. It supports giant dipole states where the electrons' center of 
mass is aligned along the electric field, with all electrons in a symmetric circular configuration 
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in the orthogonal plane. These states have been investigated for an arbitrary number of 
electrons in a normal-mode analysis, and numerically for the case of two-electron giant 
dipole states employing wave-packet propagation. 

The normal-mode analysis of the A^-electron atoms leads to a quadratic eigenvalue prob- 
lem. Its eigenmodes (as well as the eigenvectors) have been studied depending on the field 
strengths B, E as well as the electron number N . A subsequent classification according to 
their characteristics indicated that the ?>N modes were grouped into N cyclotron modes and 
N Coulomb modes, corresponding mostly to the motion induced by the magnetic field and 
the Coulomb potential, respectively. Moreover, we found a CM mode reflecting the center-of- 
mass gyration in the magnetic field, as well as a zero mode referring to a rotational invariance 
of the extremal configuration, also studied in detail. The residual modes were termed decay 
modes on account of their peculiar behavior suggesting instability for atoms with three or 
more electrons. Their absence for the system N — 2 provided a strong indication of local 
stability, with the single exception of the so-called horizontal configuration. 

In the second part of the present work, an ab-initio simulation of the six-dimensional 
two-electron system was performed. Since this type of investigation can be viewed as the 
state of the art of what is numerically feasible, we resorted to the Multi-Configuration 
Time-Dependent Hartree (mctdh) method, a wave-packet dynamics tool known for its 
unique efficiency in higher dimensions. Both the stability of the six degrees of freedom and 
some spectral properties have been examined. To establish the link to the normal-mode 
analysis, the influence of the rotational freedom of the circular conflguration has also been 
investigated. We flnd that the stability of the system strongly increases for larger electric 
flelds. For a fleld strength twice the critical value, some modes experience a decay on the 
time scale of 10"^ — 10^ ps; for the tenfold critical value, there is no instablity on the scale we 
considered, i.e., 10^ ps. These states proved rather robust against perturbations, simulated 
by displacements parallel to the magnetic fleld, where the motion is not gyrationally 
stabilized. As anticipated, the stability of the giant dipole states turned out to apply only 
outside the singular 'horizontal conflguration' of the extremum. 

While the investigation into the local aspects of stability can be regarded as somewhat 
completed, a rigorous numerical analysis is to date limited with an eye toward time and 
computational effort. First and foremost, the mctdh calculations set up so far could in 
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principle be modified to the case of > 3 electrons. However, this appears to be a massive 
numerical challenge, since the mctdh method is slowed down significantly the more degrees 
of freedom are involved. One approach to circumvent this would be an adiabatic separation 
similar to [? ], amounting to an average over the rapid cyclotron motion. This would allow 
for a significantly faster integration. For future investigations, it may also be beneficial 
for short-time propagations to implement a discrete- variable representation of the Landau 
orbitals in mctdh. 

Apart from these questions, concerned with a wider range of application and efficiency, 
the task arises to make the analysis more quantitative. A systematic resonance calculation 
yielding both the energy levels and decay widths of the giant-dipole states is desirable, but 
has so far proven to be a severe challenge computationally due to the type of instability. It 
might also yield a deeper understanding of the resonance wave functions. Ultimately, one 
could even go so far as to simulate the preparation of giant dipole states, migrating from the 
Coulomb well to the decentered extremum. In analogy to the scheme suggested for single 
electrons [? ], this may serve as a bridge to the experimental verification of giant dipole 
resonances. 
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